The Equation of State for Dense QCD and Quark Stars 
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We calculate the equation of state for degenerate quark matter to leading order in hard-dense- 
loop (HDL) perturbation theory. We solve the Tolman-Oppenheimer-Volkov equations to obtain 
the mass-radius relation for dense quark stars. Both the perturbative QCD and the HDL equations 
of state have a large variation with respect to the renormalization scale for fi < 1 GeV which leads 
to large theoretical uncertainties in the quark star mass-radius relation. 
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I. INTRODUCTION 

An understanding of the behavior of quantum chromo- 
dynamics (QCD) at high density is crucial to describing 
the physics of compact stars. This is due to the fact 
that the nuclear matter within these compact objects 
may be sufficiently dense to undergo a phase transition 
to a deconfmed or "quark-matter" phase. Of particular 
interest is the possibility that stars which have a signifi- 
cant quark-matter component could have a dramatically 
different mass-radius relationship than normal neutron 
stars jjj. However, in order to make definitive statements 
about the mass-radius relationship for such objects one 
needs a reliable calculation of the equation of state of 
high-density QCD. 

In recent years high-density QCD has received con- 
siderable attention due to the possible breaking of color 
gauge symmetry which gives rise to color superconduc- 
tivity. As a consequence of this, QCD has a very compli- 
cated phase diagram which depends on the number and 
masses of the dynamical quarks ||. It would therefore 
seem that the effect of the superconducting phase on the 
equation of state would need to be taken into account; 
however, despite the fact that the gaps for color super- 
conductivity are inherently non-perturbative, their effect 
on the equation of state for high-density QCD is expected 
to be small ||. Therefore, it is possible to use equations 
of state that have been derived neglecting them. 

The canonical choice for the equation of state for 
quark-matter has been to apply non-ideal bag models 
with various values for the bag constant M. However, 
others have applied quasiparticle models raj or directly 
applied the QCD weak-coupling expansion to first or sec- 
ond order in the strong coupling constant a s Q. The 
possibility of applying the weak-coupling expansion is an 
enticing option since this expansion can be unambigu- 
ously derived from first principles. However, since the 
strong coupling constant is expected to be on the order 
of one in the phenomenologically relevant density range, 



the question of the convergence of the weak-coupling ex- 
pansion becomes a very important one. 

At high temperatures, for example, the weak-coupling 
expansion for the QCD equation of state has been 
carried out to order a b J 2 Unfortunately, the 

finite-temperature weak-coupling expansion converges 
very slowly. In order to improve the convergence of 
the weak-coupling expansion calculational frameworks 
based on hard-thermal-loop resummation have been pro- 
posed P-p^l- These approaches attempt to describe 
finite-temperature QCD in terms of weakly-interacting 
massive quasiparticles by reorganizing the perturbative 
expansion around a state which includes hard-thermal- 
loop quasiparticles at lowest order. Detailed studies us- 
ing these techniques have shown that the convergence 
of the reorganized perturbative series appears to be 
much better than standard resummed perturbation the- 
ory HQ. 

Motivated by the success of the hard-thermal-loop re- 
organization of perturbation theory we apply the equa- 
tion of state for dense QCD at zero temperature obtained 
from hard-dense-loop perturbation theory (HDLpt) to 
quark stars. The goal of doing this is to see if this tech- 
nique can reduce the scale-dependence of the final results 
and improve the convergence of the successive approxi- 
mations to the finite-density QCD equation of state. In 
this paper we calculate the equation of state of quark- 
matter to leading order in HDLpt and use the result to 
calculate the mass-radius relation for non-rotating quark 
stars. We also make comparisons with the QCD weak- 
coupling expansion and discuss the large theoretical un- 
certainties in the weak-coupling and HDLpt results. 

The paper is organized as follows. In section II, we 
discuss the weak-coupling expansion equation of state. 
In section III, we list the finite-temperature and den- 
sity expressions for the leading order HTLpt /HDLpt free 
energy along with the HTL/HDL quark and gluon self- 
energies. In section IV, we derive the leading order 
HDLpt equation of state at zero temperature and finite 
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density. In section V, we use the resulting HDLpt equa- 
tion of state to determine the mass-radius relationship 
of a non-rotating quark star. Finally, we summarize in 
section VI. Necessary integrals are tabulated in the Ap- 
pendix. 



II. WEAK-COUPLING EXPANSION 

The zero-temperature, finite-density weak-coupling 
expansion for the free energy of an SU(N C ) gauge 
theory with Nf massless quarks has been calculated 
through order a 2 s by Freedman and McLerran [ fl5| |, and 
by Baluni [ fl6| using the momentum-space subtraction 
scheme. In the modified minimal subtraction or MS 
scheme the free energy for iV c = 3 is 
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where A = A -0.536N f +N f \ogN f , A = 10.734±0.13, 
A is the renormalization scale, and /i is the quark chem- 
ical potential 0|| . 

For the scale dependence of a s we use the three-loop 
running 
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where f3 = 11 - 2N f /3, Pi 



51 - 19A//3, and /3 2 = 
2857 -5033A//9 + 325 Nf/27 ]lj. As the boundary con- 
dition for the integration we fix a s = 0.1181 and Nf = 5 
at A — Mz = 91.1182 GeV and decrease Nf by one as 
the bottom (M b = 4.0 GeV) and charm (M c = 1.15 GeV) 
quark thresholds are crossed. Below the charm mass we 
continue the integration with Nf = 3. Using this method 
we obtain a s (l GeV) = 0.4586. Note that this method 
is only approximate since within the MS renormalization 
scheme there are non-trivial matching conditions which 
need to be imposed as each quark threshold is crossed; 
however, the corrections are very small so we have ig- 
nored them and simply required continuity of a s | p8[ . 

In Fig. [l], we show the pressure (V = —J 7 ) of a de- 
generate quark-gluon plasma truncated at order a s and 
a 2 s . The bands shown are obtained by varying the renor- 
malization scale A by a factor of two around the central 
value of A = 2/j. As can be seen from this figure, there is 
a large theoretical uncertainty in the pressure resulting 
from the choice of the scale, particularly for [i < 1 GeV 
which is the phenomenologically relevant range. For ex- 
ample, we see that at NLO the quark chemical potential 
at which the pressure vanishes varies from 200 MeV to 
650 MeV. 
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FIG. 1. Perturbative result for the pressure of a degener- 
ate quark-gluon plasma as a function of the quark chemical 
potential /x truncated at order a s and o? a . Bands correspond 
to variation of the renormalization scale \i < A < 4/x. 

Note that the perturbative expansion of the thermody- 
namical potential is an expansion in a s with coefficients 
that are polynomial in log(a s ) and that the logarithms of 
a s are due to the effect of plasmons. To properly assess 
the convergence of the weak coupling expansion of the 
thermodynamics potential at T = and fi 0, we need 
to determine the order a 3 contribution. While this would 
be a very complicated calculation, it can be obtained 
since the entire power series in a s can be calculated using 
diagrammatic methods (n]]2(]] . This is in contrast to the 
high-temperature case where non-perturbative methods 
are required at order a 3 s due to infrared divergences asso- 
ciated with the screening of static magnetic gluons pi] . 



III. HDL FREE ENERGY 

The weak-coupling expansion is an expansion about 
an ideal gas of massless particles and the lack of conver- 
gence suggests a reorganization of the perturbative series 
to improve the convergence and reduce the scale depen- 
dence. One possibility is to use HDLpt which is the ana- 
log of hard-thermal-loop perturbation theory (HTLpt) 
in the case of large chemical potentials. In HDLpt, 
one uses effective propagators and vertices that include 
plasma effects such as propagation of massive quasipar- 
ticles, screening, and Landau-damping. Thus the expan- 
sion point of HDLpt is that of ideal gas of massive parti- 
cles. In the case of zero chemical potential and high tem- 
perature, this way of reorganizing the perturbative ex- 
pansion has dramatically improved the convergence and 
reduced the renormalization scale dependence [p|4l2[ . 
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At one-loop, HDLpt approximates the high-density 
phase of QCD by a gas of non-interacting massive quasi- 
particles. The one-loop HDL free energy for an SU(N C ) 
gauge theory with Nf massless quarks is 



T=(N 2 -l)[(d-l)T T + T L ] 

+N c NfT q + A £ 



(3) 



where d = 3 — 2e is the number of spatial dimensions and 
Tt and Tl are the contributions to the free energy from 
the transverse and longitudinal gluons, respectively ]l0| |. 
T q is the contribution to the free energy from each color 
and flavor of the quarks, and Ao£o is the leading-order 
vacuum-energy counterterm. The contributions from the 
transverse and longitudinal gluons are 



?t = ^io g [p 2 + n T (p)] , 
^ = ^io g [p 2 -n i (p)] , 
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where Ht(P) and Hl(P) are the transverse and longitu- 
dinal HDL gluon self-energies. The quark contribution 



T q is 



logdet [f-S(P)] 



(5) 



where S(P) is the quark self-energy. The quark contri- 
bution can be rewritten as [BTiEl 
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The functions II T (P), U L (P), A (P), and A S (P) are de- 
fined by 
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where m g and m q are the gluon and quark mass param- 
eters respectively, and the function 7p is defined by 



T P = w(e) / dc (1 - c 2 Y 
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where the function w(e) is 
w(e) 
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IV. ZERO TEMPERATURE AND FINITE 
CHEMICAL POTENTIAL 

In this section, we calculate the zero-temperature limit 
of the one- loop HDL free energy Eq. (||). The zero- 
temperature limit of the gluon contribution was calcu- 
lated in Ref . jl(| , while the the zero-temperature limit of 
the quark contribution was calculated in Ref. p2| ] using 
three-dimensional expressions for the self-energies. Since 
we are using dimensional regularization in our calcula- 
tions, we use d-dimensional expressions for them. To 
leading order in HDL perturbation theory, the difference 
can be absorbed in a redefinition of the renormalization 
scale A. For higher order calculations, it essential to use 
the d-dimensional expressions for the quark and gluon 
self-energies. 



A. Gluon contribution 

Our original expression for the gluon contribution 

T g = {d- 1)T T + Fl (13) 

to the free energy involves a sum over the discrete Mat- 
subara frequencies u> n — 2irnT. As T — > 0, the sum ap- 
proaches an integral over the continuous energy to. The 
only scale in the resulting integral is m g and T g is pro- 
portional to rrig on dimensional grounds. 

In the zero-temperature limit, the transverse free en- 
ergy Tt becomes 
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Since Hp is a function of the combination iv/p only, it 
is convenient to rescale the energy uj — > pto. Integrating 
over the angles of p and using the fact that the integrand 
is an even function of u>, the integral reduces to 
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where fid = 2-K d / 2 /T(d/2) is the angular integral. The 
dimensionally regulated integral over p can be evaluated 
analytically using Eq. (A. 9) giving 
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Expanding around d = 3 and evaluating the resulting 
integral over u> numerically, we obtain 
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T T = - 
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The longitudinal free energy Tl can be evaluated in the 
same manner and reads 
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0.320878 



(18) 



The gluon contribution J-" g to the free energy is obtained 
by inserting Eqs. (0) and (|l|) into Eq. Q 
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1.24546 



(19) 



The pole in Eq ( pL 9|) agrees with the one found in Ref. [|10| , 
but the finite term differs since d was set equal to 3 in 
the expression for Tp 



B. Quark contribution 



The total free energy is given by the sum of 
Eqs. ©, ©, and (H): 
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Our leading order result for the thermodynamic func- 
tions depends on the gluon and quark mass parameters, 
and the renormalization scale A. These parameters are 
completely arbitrary in the sense that the dependence 
on them will be systematically subtracted out at higher 
orders. If higher order calculations were available, the 
masses could be determined by a variational principle 
giving rise to a self-consistent gap equation |]Jl^]. In 
a one-loop calculation, we have little other choice than 
taking the weak-coupling expansion results for the gluon 
and quark masses. The gluon and quark masses are p3| 
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The quark contribution Eq. Q can be expanded in a 
power series of m 2 //i 2 . To second order in to 2 //! 2 , we 
obtain 
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Using the results for the zero-temperature limit of the 
different sum-integrals listed in the Appendix, Eq. (|2T 
reduces to 
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We note that the quark contribution to the free energy 
is ultraviolet finite. The coefficient of the to 4 , term is dif- 
ferent from the one found in Refs. |h],|22|], since d was set 
equal to 3 in the expression for Tp. 



C. One- loop free energy 

The gluon contribution to the free energy Eq. (|n]) is ul- 
traviolet divergent, while the quark contribution Eq. ( |2p| ) 
is finite. The ultraviolet divergence is cancelled by the 
counterterm that was determined in Ref. fic| l 
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In the remainder of the paper we specialize to case 
N c — Nf — 3. The free energy (|23] ) then reduces to 
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Comparing the weak-coupling expansion Eq. ([I]) with 
the leading order HDL result Eq. (^3|), we note that 
the order-a s term is over-included by a factor of two 
but the a 2 log a s term is included exactly. A next-to- 
lcading order calculation in HDL perturbation theory 
would agree with the weak-coupling expansion at orders 
a s and a 2 log a s if we identify the gluon and quark mass 
parameters with their weak-coupling expressions. 

In Fig. ||, we show the leading-order HDL result for 
the free energy T normalized to that of an ideal gas. For 
comparison, we also show the next-to-leading order pre- 
diction from the QCD weak-coupling expansion as a grey 
band. The band is obtained by varying the renormaliza- 
tion scale A by a factor of two around the central value 
A = 2/.1. In this figure, we see that the HDLpt results 
also have a large variation with respect to the renormal- 
ization scale. Additionally, we see that the results for 
both A = 4/i and A = 2/i are both unphysical in that 
they predict negative quark number densities. If we add 
a requirement that the quark number density be positive 
we find that this requires A < I.6/1. We will use this as 
the lower bound for A in the plots of the mass-radius 
relation in the next section. 
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FIG. 2. Leading-order HDLpt result for the pressure of a 
degenerate quark-gluon plasma as a function of chemical po- 
tential (i. The NLO weak-coupling expansion result is shown 
as a grey band. Band corresponds to variation of the renor- 
malization scale fi < A < 4/x. 

Note that the requirement that A^1.6/i may have 
some physical basis since the scale of the coupling con- 
stant should be related to the average momentum ex- 
change of two quarks on the Fermi surface. At zero tem- 
perature the largest momentum exchange possible is 2fi 
and the smallest momentum exchange is of the order of 
the superconducting gap (p. Therefore, the scale for the 
coupling constant should be in the range <, A < 2jU so 
that the choice of A ~ 1.6/1 is not unreasonable. 



V. MASS-RADIUS RELATIONSHIP 

The mass-radius relationship for a non-rotating spher- 
ically symmetric star is obtained by solving the Tolman- 
Oppenheimer-Volkov (TOV) equations |24|] for the mass 
M and the pressure (V — ~T) as a function of the radial 
distance from the center: 



= Airr 2 E(r) 

dr 



dp 
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2GM(rY~ 1 



(27) 



(28) 



where G is Newton's constant, c is the speed of light, 
£ = £/c 2 , andV = V/c 2 . 

In this work we will ignore the presence of the nuclear 
phase of matter which is expected to undergo a first-order 
phase transition to the quark-matter phase. A more de- 
tailed study would include the effects of the nuclear phase 



on the mass-radius relationship; however, our goal here 
is only to show that both standard perturbation theory 
and HDLpt have large theoretical uncertainties related to 
the renormalization scale dependence. The most plausi- 
ble scenario is that there will not be "naked" quark stars, 
but instead there will be neutron stars with a very com- 
pact quark-matter core and a thick outer layer of normal 
nuclear matter. 
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FIG. 3. Mass-radius relation for a quark star with 
A//i = 1.6 and = 1. The weak-coupling results for the 
same choice of renormalization scales are shown as dashed 
lines. Mq — 1.989 x 10 30 kg is the mass of our sun. 

In Fig. ^, we show the mass-radius relationship ob- 
tained by solving the TOV equations numerically for 
A//x = 1.6 and A/// = 1. For comparison, we also 
show the QCD weak-coupling expansion results for the 
same choice of renormalization scale as dashed lines. As 
can be seen from this figure there is a large variation 
in the mass-radius relationship as the renormalization 
scale is varied over even this rather limited range of 
fj, < A < 1.6 fi. Using this range, we find that using 
the HDLpt equation of state ( p6|) that i? max ~ 3.4 — 10.9 
km and M max ~ 0.6 — 2.12 Mq. With this same range 
we find that using the perturbative equation of state (0) 



that i? n 



2.4 - 5.6 km and M n 



VI. DISCUSSION 



0.42 - 0.95 M, 
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In this paper, we have calculated the free energy of cold 
dense quark matter to leading order in HDL perturbation 
theory (HDLpt). The predictions of HDLpt depend on 
a renormalization scale A that arises both from running 
of the coupling constant and from the renormalization of 
the additional ultraviolet divergences that are introduced 
by the HDLpt reorganization of perturbation theory. It is 
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possible to separate these two effects by introducing two 
renormalization scales A3 and A4 as done in Ref. Jl0| , 
and these scales are associated with the soft scale and 
hard scale, respectively. For simplicity we have chosen 
not to distinguish between the two and have simply set 
A 3 = A 4 . 

We then used the HDLpt equation of state as input to 
the TOV equations in order to determine the mass-radius 
relationship of quark stars. We find that the large scale 
dependence of the HDLpt and weak-coupling expansion 
equations of state lead to large theoretical uncertainties 
in the quark star mass-radius relationship. In the case of 
the weak-coupling expansion, the expansion only seems 
to be under control for fj, > 5 GeV. For HDLpt, the scale 
dependence is larger for all values of fi and without a 
next-to-leading order HDLpt calculation it is not pos- 
sible to draw conclusions about the convergence of the 
series. In addition, the choices A = 4/x and A = 2/i lead 
to rather unphysical predictions as can be seen in Fig. ||. 
In order to eliminate these we were forced to further re- 
strict the range of renormalization scales considered to 
/i < A < 1.6 fi. The failure of both HDLpt and the weak- 
coupling expansion to reliably describe the finite-density 
QCD equation of state for /i between 300 MeV and 1 
GeV is troubling since this is the range which is impor- 
tant for determining the mass-radius relationship for a 
quark star. 

As mentioned in Section [□], it possible that a compu- 
tation of the order contribution to the finite-density 
QCD equation of state could remove some of the theo- 
retical uncertainties resulting from the use of the weak- 
coupling expansion result. Despite the fact that this 
would be a rather difficult task it seems that this cal- 
culation is required in order to draw more firm conclu- 
sions about the QCD equation of state. However, even 
if this calculation were available, the presence of a non- 
perturbative contribution from a color-superconducting 
phase of QCD in this range of quark chemical poten- 
tial introduces additional theoretical uncertainties. Per- 
turbative results extended down to this range of quark 
chemical potential give gaps on the order of cj) ~ 30-100 
MeV. Since the gap gives a relative contribution of the 
order of <fi 2 / ji 2 this could translate into a relative modi- 
fication of the equation of state between 1% and 10%. 

A very challenging problem would be to calculate the 
next-to-leading order correction to the free energy in 
HDLpt. If the next-to- leading order correction turns out 
to be small for relevant values of the chemical potential, 
the results obtained in the present work can be trusted. 
However, it would seem more prudent to compute the or- 
der a' contribution in the weak-coupling expansion since 
the finite-density perturbation series does not seem to 
suffer from the same problems (oscillation and lack of 
convergence) as the finite-temperature perturbation se- 
ries. 
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APPENDIX A: 

In the imaginary-time formalism for thermal field the- 
ory, the 4- momentum P = (Po,p) is Euclidean with 
P 2 = P 2 + p 2 . The Euclidean energy po has discrete 
values: Pq = 2nirT for bosons and Pq = (2n + 1)ttT + /i 
for fermions, where n is an integer and /x is the chemi- 
cal potential. Loop diagrams involve sums over Pq and 
integrals over p. With dimensional regularization, the 
integral is generalized to d = 3 — 2e spatial dimensions. 
We define the dimensionally regularized sum-integral by 
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where 3 — 2e is the dimension of space and A is an arbi- 
trary momentum scale. The factor (e 7 /47r) e is introduced 
so that, after minimal subtraction of the poles in e due 
to ultraviolet divergences, A coincides with the renormal- 
ization scale of the MS renormalization scheme. 



1. Simple one-loop sum-integrals 

The simple fermionic sum-integrals required in our cal- 
culations are 
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These sum-integrals can be calculated by standard con- 
tour methods. 
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2. One-loop HDL sum-integrals 

The one-loop fermionic sum-integrals involving the 
HDL function Tp are 
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These sum-integrals can be calculated using the methods 
developed in Rcf. (Til. 



3. Integrals 

In order to calculate the zero-temperature limit of T g , 
we need the following integral 

f°° p ( a+1 \ y I l — a \ 

/ dpp a log{p 2 +m 2 ) = V 2 ; V 2 ' m a+1 . (A.9) 
7o ' " + 1 
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